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ABSTRACT 

We examine the tidal disruption event scenario to explain Sw 1644+57, a powerful and persistent X-ray 
source which suddenly became active as GRB 110328A. The precise localization at the center of a z = 0.35 
galaxy argue for activity of the central engine as the underlying cause. We look at the suggestion by Bloom 
et al of the possibility of a tidal disruption event (TDE). We argue that Sw 1644+57 cannot be explained 
by the traditional TDE model in which the periastron distance is close to the tidal disruption radius — three 
independent lines of argument indicate the orbit must be deeply plunging or else the powerful jet we are 
observing could not be produced. These arguments stem from (i) comparing the early X-ray light curve to 
the expected theoretical fallback rate, (ii) looking at the time of transition to disk-dominated decay, and (iii) 
considering the TDE rate. Due to the extreme excess in the tidal force above that which would be required 
minimally to disrupt the star in a deeply plunging orbit at periastron, we suggest this scenario might be referred 
to more descriptively as a TOE (tidal obliteration event) rather than a TDE. 

Subject headings: accretion, accretion disks - black hole physics - galaxies: active - galaxies: nuclei 



1. INTRODUCTION 

Given the massive black holes (MBHs) known to reside at 
the centers of most galaxies (Magorrian et al. 1998), an in- 
evitable occurrence will be the occasional errant stellar wan- 
dering into the inner galactic environment, where the prospect 
for a tidal disruption and subsequent swallowing of some frac- 
tion of the ill-fated star may ensue (Rees 1988, 1990). A 
complete study of the tidal disruption process was carried 
out recently by Strubbe & Quataert (2009, hereafter SQ09) 
and Lodato & Rossi (2011, hereafter LRU). For MBHs with 
greater than 10^ Mq, the tidal disruption radius lies inside the 
event horizon for a typical main sequence star, so the star 
would be swallowed whole. The rate of such events is cal- 
culated to be ~10~^ — 10"'^ yr~^ galaxy"^, depending on 
galaxy type (Wang & Merritt 2004). Wang & Merritt de- 
rive a rate ^10^"'' yr^^ Mpc^'^, which, given a local galaxy 
space density of ^0.01 Mpc^^ (Driver et al. 2005), yields 
an effective rate per galaxy of about 10~^ yr^^. This is ^10 
times higher than previous estimates (e.g., Cohn & Kulsrud 
1978). The best candidates to date are presented in van Velzen 
et al. (2010) — two strong examples of optical flares from 
tidal disruption events in archival SDSS data, from which 
the authors determine a rate of tidal disruption events (TDEs) 
^3^2 X 10~^ yr~^ galaxy"^, about a factor of 30 below the 
Wang & Merritt theoretical estimate. There were also earlier 
claims of TDEs (e.g., Komossa et al. 2004, 2008; Gezari et 
al. 2008, 2009). 

After the star comes undone, roughly half of the remnants 
remains on bound, highly eccentric elliptical orbits and ac- 
cretes; the other half lies on unbound, hyperbolic orbits and is 
ejected (Lacy, Townes, & Hollenbach 1982, hereafter LTH82, 
Rees 1988, 1990). Theoretical research on accretion of tidal 
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debris from tidal disruptions has traditionally been divided 
into two phases: (1) the immediate unbinding of the star 
and accretion of streamers of gas, for which estimates a fall- 
back rate M ^ t'^/^ (Rees 1988; Evans & Kochanek 1989; 
Lodato, King, & Pringle 2009, SQ09, Ramirez-Ruiz & Ross- 
wog 2009, Guillochon et al. 2009, LRU), and (2) the sub- 
sequent disk accretion, for which the freely expanding outer 
disk enforces M oc t~^^/^^ for a thin disk (Cannizzo, Lee, & 
Goodman 1990, hereafter CLG90, Ulmer 1999, SQ09). and 
M a t-^/^ for a thick disk (Cannizzo & Gehrels 2009, here- 
after CG09). The longer term evolution, which can be cast in 
terms of accretion disk physics, predicts the tidal disruption 
flare to be a strong UV source, ~10*^ erg s~^, which a du- 
ration of months to years. These expectations are nominally 
satisfied in the handful of secure candidates (van Velzen et al. 
2010). 

Sw 1644+57 was triggered by the Swiftl^P^ as a low SNR 
image trigger requiring a 2000s integration (Cummings et al. 
201 1 ; see Burrows et al. 2011 for full details of the Swift anal- 
ysis). Ground-based NIR and radio observations provided the 
first clear evidence connecting the X-ray/gamma-ray source 
with the galaxy nucleus (Berger et al. 201 1). Precise localiza- 
tions by HST (Fruchter et al. 201 1) and Chandra (Levan et al. 
201 la) confirmed the transient to he at the center of its galaxy. 
Sw 1644+57 is at a redshift z = 0.3534 ± 0.0002 (Levan et 
al. 201 lab). The inferred R-band magnitude at that distance 
implies Mr ~ —18 (Fruchter et al. 2011), which implies a 
~1O^M0 black hole at the galactic center, if the galaxy fol- 
lows the Magorrian relation (Magorrian et al. 1998). Burrows 
et al. (2011) derive a central BH mass 2 x lO^'Mp based on 
the Magorrian relation, and a lower Umit~7 x 10" Mq based 
on X-ray variability. Quataert & Kasen (201 1) point out how- 
ever that the signal-to-noise of the XRT data is not sufficient 
to constrain significant variability on < 10 s timescales, there- 
fore the central BH mass could be much lower. Similarly, the 
host galaxy brightness places only an upper limit on the BH 
mass from the Magorrian relation, since no "bulge" is actu- 
ally resolved. In addition. Miller & Giiltekin (2011) use em- 
pirical relations relating central black hole mass to radio and 
X-ray luminosity for Swl644+57 to argue for a central mass 
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log(MMBH/-W0) = 5.9 ± 1.1. In view of these considera- 
tions and the uncertainty associated with applying the Magor- 
rian relation, in this work we seriously entertain the possibiUty 
that Mmbh may he between ~ 10*^ Mq and IO'^Mq. 

The positional coincidence gives strong credence to the no- 
tion of an event associated with the putative MBH at the 
galaxy center, quite probably a tidal disruption event (Bloom 
et al. 201 lab, Shao et al. 201 1), and yet the contrast between 
theoretical expectation (UV source at 10"''' erg s^^) and obser- 
vation (gamma-ray/X-ray source at 10^'' — 10^^ erg s~^) lead 
to the unavoidable conclusion that we are viewing a strongly 
beamed event. In fact, the jet that we must be observing is 
completely divorced from the physical properties of the fall- 
back disk, other than relying on accretion as a power supply. 

In this work we estimate the rate of accretion within a tran- 
sient disk formed from the tidal debris of a disrupted star, and 
estimate a temperature in the inner disk. We also calculate the 
expected rates, both for SwiftfBAT and aU-sky X-ray monitors 
with ~mCrab sensitivity. 

2. GENERAL FRAMEWORK 
2.1. Accretion Properties 

Given the existence of a compact mass Mmbh and the 
sudden introduction of a much smaller mass of stellar de- 
bris AM at ~ 1 — IO-R5, where the Schwarzschild radius 
Rs = 2G'Mmbh/c^, what will be the subsequent rate of ac- 
cretion onto the MBH as a function of time? 

In simple physical terms, the disruption of the star occurs 
when the star comes closer to the MBH than a tidal disrup- 
tion radius Rt which is determined by demanding that the 
mean density of the volume of space enclosed by Rt, i.e., 
AfMBH/(47ri?y/3), be equal to the density of the star. This 
gives 

Rt ~ i?*(MMBH/M,)i/3 = 1.50 X 10" cmm-'^\,ml^\ 

(1) 

where = M^,/Mq, = R^,/Rq, and = 

Mmbh/1O^M0 (Rees 1988, 1990). Expressing Rt in terms 
of Rs = 2.95 X 10^^ cm mb7 gives 



Rt 



r 1 -1/3 -2/3 

5.1 m« r*mj7 , 



(2) 



Setting this equal to unity and solving for mbr, we see that for 
mb7 > 10, the tidal radius will lie inside the Schwarzschild 
radius, and a main sequence star will be swallowed whole. For 
a maximal Kerr BH in which Rs = GMmbh/c^, there would 
still be another factor 2^/^ increase in A/mbh before this limit 
is reached. An equivalent way of visualizing this is in terms 
of the mean density enclosed by Rs- For a nonspinning (i.e., 
Schwarzschild) BH, the mean density out to Rs is pmbh = 



185 g cm 
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For comparison Pq = 1.4 g cm Stellar 



disruption occurs if pmbr > P* and Rs < Rp < Rt, where 
iip is the periastron distance of the star. 

We define ^ = Rp/Rt- The requirement that the TDE 
occur outside the event horizon enforces 



^ > xRs/Rt « 0.2 rnj\:^m%''x, 



(3) 



where x = 1 for a Schwarzschild BH and x = 0.5 for a Kerr 
BH. 

We now look at the long term evolution of the resultant 
mass accretion rate onto the MBH, both in terms of stellar 
debris fallback and accretion disk. 



2.1.1. Debris from Stellar Fall-Back 

Earlier studies give an expression, derived in the Newtonian 
limit, for the time scale of return of the most bound stellar 
material to periastron (LTH82, Rees 1988, cf. Equation [4] of 

LRU) 

TT (RpV''' 
ifallback — -^Yp. \ j 

where the dynamical time at periastron tp = 
(GMmbh/-Rp)~^/^. Scaling Rp in terms of Rr gives 

ifaiiback = 130 d ml^^^^m^^r^^'^. (5) 
The condition on ^ from Equation (3) gives the lower limit 

tfaiiback > 0.98 d x^ml^^r-^^^ ■ (6) 

The accretion rate due to fallback is then (e.g., LTH82, Rees 
1988, SQ09,LR11) 



Mi, 



1 



"allback — 



3 tfallback V ^fallback 



t 



-5/3 



(7) 



so that Mfaiiback oc m^^^^^. The total fallback mass 
AM{ti,t2) = Mfaiibackrfi 0.5M* in the limit ti 

^fallback and ^2 -> 00. Note that Affaiiback = for t < 
^fallback, and we can thus identify t = tfaiiback as the time 
of the earliest significant BAT activity, ^ 2 d before the GRB 
110328A trigger. 

We can estimate the basic dynamical parameters of the 
tidal disruption event in the following way. Superposed on 
the long-term X-ray light curve of Sw 1644H-57, Levan et al. 
(2011b) and Bloom et al. (2011b) show a t'^'^ decay law. 
This suggests that the X-ray luminosity tracks the rate of re- 
turn of the fallback material to pericenter. In this work how- 
ever, we argue that only the early light curve, i < 6 d, is repre- 
sentative of the stellar debris fallback; the subsequent evolu- 
tion is dominated by the accretion disk. As mentioned above, 
the very high X-ray luminosity points to a strongly beamed 
source, with a beaming factor Cb ^ 0.01, so that the intrinsic 
jet luminosity ij is related to the observed Lx by 



(8) 



Let us assume that the jet luminosity is a fraction / of the 
kinetic energy carried by the outflowing material: 



Li 



(9) 



where we have also assumed that the outflow velocity ~ c. 
The outflow rate itself is a fraction ej of the fallback rate 
Mfaiiback^ where, for a mildly super-Eddington accretion flow, 
we expect ej - 0.1 (SQ09, LRU, Dotan & Shaviv 201 1). We 
do not wish to restrict ourselves to baryonic jets; this same 
analysis could be apphed to a Poynting flux jet, with / de- 
noting a mass-energy equivalent flux. We can thus obtain an 
estimate of the fallback rate Mfaiiback as a function of the ob- 
served X-ray luminosity: 



Mfaiiback (i) = 



hi c2 ■ 



(10) 



Using Equations (5) and (7), we can thus obtain an estimate 
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of the fallback time and of the penetration factor ^: 



fallback 



(11) 



t 



130d 



1/3 



1/2 



1/3 -1/2 -1/6 



(12) 

Note that ^ is very weakly dependent on both the stellar mass 
and on the black hole mass £_ cx (m*TOh7)-i/6 (Krolik & Pi- 
ran 2011). We acknowledge the simplicity of the assumption 
that / and and ej are independent of the Eddington ratio of the 
accretion rate. There is observational evidence that AGN de- 
crease in radio loudness with decreasing Eddington ratio (e.g., 
Sikora et al. 2007). This may be the result of a fundamental 
change in the underlying accretion disk structure, e.g., from a 
slim disk to a thin disk. Based on the total Sw1644h-57 light 
curve to date, which has exhibited a dynamic variation in av- 
eraged X-ray flux of <2 orders of magnitude, our assumption 

/ = e'j = e'b = may be at least crudely justified. In other 
words, although the accretion physics will likely change with 
Eddington ratio (i.e., slim-to-thin disk), this transition occurs 
at a low Eddington ratio ^0.01, which is not applicable to 
date for Sw1644h-57. 

Several works have considered the observational conse- 
quences of ultra-relativistic jets from TDEs prior to Sw 
1644H-57 (Farrar & Gruzinov 2009, Giannios & Metzger 
2011). Giannios & Metzger find, on the assumption of a 
baryonic jet, that emission from the forward shock may be de- 
tectable for several years. However, motivated by results from 
recent groups studying 3D GRMHD jets (Hawley & Krolik 
2006, McKinney & Blandford 2009) which find a baryonic 
zone-of-exclusion in the jet beam, leading to primarily Poynt- 
ing flux jets, we do not restrict our focus to baryonic jets, and 
indeed we view Poynting flux jets to be more likely. On the 
other hand, the jet luminosity can be dominated by Poynting 
flux even if most of its rest mass is in baryons. If the jet re- 
mains Poynting dominated out to the distance where it begins 
to interact with the ISM, a forward shock with concomitant 
non-thermal emission can still be produced (e.g., Mimica, Gi- 
annios, & Aloy 2009). A highly magnetized jet transfers en- 
ergy efficiently to the environment, and the Lorentz factor of 
the forward shock is related to the magnetization of the jet, 
not its Lorentz factor, which could be much less. 

At t = 2.5 d from the BAT trigger the X-ray luminosity 
crosses 10^^ erg s^^. Inserting these values into Equation 
(11), after taking into account the offset time 2 d to be added 
to the 2.5 d, we obtain ifaiiback ^ 6480 s {(.^j < 2.5 

d. From Equation (12) we obtain an estimate of ^: 



C 0.0833 



Vo.oi/ Vo.iy 



-1/2 



0.1 



-1/2 



(13) 



strained parameters, / and e^. By requiring that ^ be larger 
than the minimum value, we get 



< 0.18x 



-2/3 -5/3 -1 



(14) 



Unless ftj/eb is much smaller than the value above, the event 
must have been deeply plunging. Note that the condition 
above is strongly dependent on the black hole mass. If the 
MBH mass were smaller, closer to 10® A/©, that would relax 
the stringency of the deeply plunging constraint. 

2.1.2. Accretion Disk 

The need for a deeply plunging orbit can also be obtained 
considering the properties of the accretion disc formed by 
the debris. The process of fallback does not directly supply 
gas into the MBH, but rather it returns material to ^Rp and 
therefore feeds an accretion disk. Given the small inferred 
^fallback valuc with rcspcct to the total length of the X-ray 
light curve, it seems probable that at early times an accre- 
tion disk is established. An evaluation of ^ -MtaiibackC^i 
using Equation (7) shows that 90% of the total available fall- 
back mass G.5Af* has already been supplied to the accretion 
disk by t2 =32ifanback- It is therefore of interest to estimate 
the properties of the accretion disk, given a resident disk mass 
-G.5A/0. 

One may estimate the long term accretion rate in Sw 
1644H-57 by looking at the properties of the accretion disk 
formed from the debris. We idealize the disk as extending 
from the event horizon to rjRp, where ry ~ 2. For the high 
rates of accretion expected we use the super-Eddington slim 
disk considered by CG09 (see their Section 2.3.2) in which 
the radial distribution of surface density S(-R) in a roughly 
steady disk would be oc i^^^/^^ thus S(r) = Y.oiR/ RpY^''^ . 
Therefore integrating the mass distribution J 2TTRdRJ^{R) 
out to rjRp, and setting it equal to the mass of the fallback 
debris AM, one obtains 



AM 



AM 



{2/3)Trr]^Rp' Svr (r,/2)^^^RT' 



(15) 



We derive an estimate for the properties of the resultant disk 
using algebraic scalings. Inverting Equation (24) of CG09 to 
obtain Mq yields 



Mn 



1.63 X 10^^ gs-^Eo™ 



1/2 1/2 
67 ^13 "-1' 



(16) 



where the disk radius ri3 — rjRp/lQ^^ cm and q;_i is the 
Shakura & Sunyaev (1973, hereafter SS73) a parameter in 
units of 0. 1 . Substituting for Eq gives 



Mo = 7.72 X 10^9 g s^^AA/i/2m^^"ri3''"a_i 



1/2 -3/2 



-1/2 



(17) 



We thus see that the numbers involved point strongly towards 
a deeply plunging evenj^with ^ <C 1. 

If one assumes 6^ = 0.01 and ej = f = 0.1, the resulting ^ 
value would be smaller than the minimum given by Equation 
(6). This allows us to put reasonable bounds on the least con- 

' Note that for a WD instead of a MS star, i.e., i?* = 5 X 10* cm and 
mi,7 = 10-^ (Krolik & Piran 2011), Equation (13) gives ^ ~ 3.1, which 
would require an adjustment of the efficiency parameters to still yield a tidal 
disruption, ^ < 1. 



where AAI1/2 is the fallback debris mass in units of 0.5Mq. 
Substituting for the disk radius r/Rp in units of lO^'^ cm, 

T r ^ -1/3 1/3 . 

ri3 — 1.5774m* r*TOj7 , gives 

Mo - 4.21 X 10^9 g s-i(r?/2)-3/2r^/'m3/Vr^/2a_i. 

(18) 

Equation (18) only gives a measure of the initial value of the 
rate of accretion within the accretion disk formed rapidly by 
stellar fallback. The longer term evolution is given by CG09 
as 



M{t) = Mo 



-4/3 



(19) 
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Fig. 1.— The SwiftfXKr light curve (Evans et al. 2007, 2009) for Sw 
1644+57, plotted on both a log-linear (first panel) and log-log scale (second 
and third panels). In the third panel we show the block-averages of the the X- 
ray flux, binned in time in 0.125 dex bins. The blue lines indicate the putative 
stellar debris fallback slope —5/3 and the red lines show the fallback decline 
rate expected due to a super-Eddington slim disk, —4/3. The data do not 
strongly support either decay law in detail, but are suggestive. 



where 



9/4 q(GMmbh)i/2 
3/2 

= 224 s ■ ' ' 



-1/2 3/2 -1 



(20) 



and tr = {GMmsyiI RtY^^^ ■ The time to is basically the 
viscous time for a thick disk at i?p. 
For a disk mass which varies as t^^ , the global viscous 

time scale tvisc(0 = | A/disk /-^^diskl — /3^^t increases with 
time. The outer edge of an "external" accretion disk, defined 
by the lack of a confining torque as from a companion in an 
interacting binary, moves outward with time. The fact that 
tvisc {t) oc t results from the spreading of the outer disk edge, 
and the viscous time is determined by the slowest time scale 
in the system, i.e., that at large radii. On the other hand, in 
comparing the time scales for mass fallback with a viscous 
time scale, the more relevant time within the disk is not that at 
the outer edge, but rather the response time of the disk at the 
radii where the fallback mass is deposited to- 

Figure 1 shows the inferred long-term X-ray luminosity for 
Sw 1644-1-57, taking into account the distance given by the 
redshift. The X-ray flux ^10^^ — 10"*^ erg s^^ is super- 
Eddington for rribj = 0.1 — 1, and also quite different spec- 



trally from theoretical expectations if we were directly ob- 
serving the disk. On the lower panels we indicate the expected 
slopes for stellar debris fallback and accretion due to the ex- 
ternal disk formed by the stellar debris. There appears to be 
a difference in the decay characteristics at t ~ 6d, which we 
interpret as the start of the disk-dominated decay. 

Figure 2 shows a comparison between the stellar debris fall- 
back rate from Equation (7) and the disk accretion rate from 
Equation (19) for ^ = 0.1 and 0.2. For mbj = 1 the tran- 
sition point for the ^ = 0.1 lines — 4.6 d is close to the 
time t ~ 6 d in the Sw 1644-1-57 light curve where there ap- 
pears to be a change in the decay characteristics. If this is 
indeed associated with the debris fallback/disk decay transi- 
tion, it provides further evidence for a small periastron value 
^ ~ 0.1. For rribj — 0.1 the putative t ~ 6 d transition cor- 
responds to a larger ^ ~ 0.18. Given the results of Miller & 
Giiltekin (201 1) the latter value may be better motivated. We 
stress, however, the uncertainties associated with 771^7 and rj, 
and therefore our inferred ^; we certainly do not claim a pre- 
cise determination of ^. On the other hand, it is noteworthy 
that the uncertainties in /, e^, and ej from the previous section 
do not enter into the (^x ) line of argument in this section, 
and yet the inferred ^ value obtained by adopting reasonable 
values for /, ef,, and ej and fitting the fallback decay to the 
early light curve are in line with that obtained using the pu- 
tative tx — 6 d value in the light curve in conjunction with 
Equations (7) and (19). 

Setting Equations (7) and (19) equal and solving for the 
time tx defining the transition point yields 

tx = 1.25 X IQi" s m:'/' r'J^ m„ a_i ivl'^f'^- 

(21) 

For our favored parameters = = a_i = 1, the transi- 
tion point in the light curve to a disk dominated decay would 
lie at 



6 



0.106 



6d 



2/9 



m 



'2/9 



(22) 



2.2. Event Rate 



Why have such events not been seen previously? As noted 
earlier, Sw 1644-1-57 was nearly missed due to its faintness in 
the ^wi/f/BAT, which is a ^mCrab instrument. Its brightness 
in the SwiftlUKY arises because of the XRT's ^/iCrab sensi- 
tivity. 

SwiftfRNT rate: The highest SNR flare from Sw 1644H-57 
seen in BAT could have been detected out to Zmax — 0.8 Q 
(Burrows et al. 2011), which encompasses a comoving vol- 
ume ^90 Gpc"^. Adopting a local galaxy density 0.01 Mpc~^, 
this gives '^lO^ galaxies in this volume. Taking a nominal 
tidal disruption rate per galaxy of 3 x 10^'^ yr^^ (van Velzen et 
al. 2010) and ^ ~ 0.1, which lowers the nominal rate by a fac- 
tor ^^■^ to 3 X 10^^ yr^^ galaxy^^, we obtain a net rate ^300 
yr^^. Adopting a beaming factor ~'10^'^ reduces the rate to 
~3 yr^^, and taking a ^wZ/f/BAT sky coverage of ^10% gives 
a further reduction to ^ 0.3 yr~^. Note that the standard TDE 
rate predicts ~'30 BAT events per year (for a beaming factor 

' The redshift ^jct = 0.35 for Sw J1644 corresponds to only ~12% of 
the estimated sensitivity volume. At present, based on a single detection, it 
is still premature to speculate on a possible cosmic evolution of TDEs. For 
the first low-redshift detection of this class of objects — a jetted TDE — the 
fact that Ztjct < 2max is most Hkely an effect of small number statistics (i.e., 
A'^ = 1). Future observations will shed light on this point. (A second jetted 
TDE candidate has recently been found with 2 = 1.2 [Cenko et al. 201 1]). 
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Fig. 2. — The theoretical accretion rates for stellar debris fallback (blue) 
and a super-Eddington slim disk (red) for the parameters given in the text. 
The two sets of lines indicate § = 0.1 and 0.2. In the upper panel we take 
mi,7 = 1 and in the lower panel mi,7 = 0.1. For rriby = 1 the transitions 
lie at t = 4.59 d ($ = 0.1, tfauback = 3.11 hr) and t = 103 d = 0.2, 
^fallback = 1-04 d), with = 0.106 corresponding to tx =6 d; for 
mi,7 = 0.1 the transitions are at t = 0.46 d (^ = 0.1, tfaiiback = 3540 
s) and t = 10.3 d (^ = 0.2, tfaiiback = 7.87 hr), with ^x = 0.177 
corresponding to t x = 6 d. 

'-^lO"^), which can be strongly excluded. This suggests that 
only a fraction of TDEs can create powerful relativistic jets 
(^10%; Burrows et al. 2011) and reinforces the notion of a 
deep plunging orbit as a way of reducing the overall TDE rate. 

Beppo-Sax rate: Beppo-Sax had a FoV ~0.12 sr (FWHM), 
a sensitivity of 2mCrab in 100 ks, a bandpass of 2 — 10 keV, 
and a lifetime 6 yr. Based on the steady state observed in XRT, 
Sw 1644+57 could have been detected up to aredshift z = 1.1 
by the Beppo-Sax/WFCs. Considering the much smaller FoV 
(0.12srvs. 1 srofBAT), the rate of events is only 0.2 in 6 yr of 
the mission, so the fact that Beppo-SAX did not see an event 
like Sw 1644+57 is consistent with the BAT rate. By using 
the brightest peaks, the event could have been observed up to 
higher redshift (z ~ 2.5), but it would have been impossible 
to identify it as a TDE because the long lasting soft emission 
would have been under detection threshold. In any event, even 
at z = 2, the increase in volume does not compensate the 
difference in FoV. 

Wide-field X-ray monitor (WFXM) rate: Future mission 
concepts with sensitive wide-field X-ray monitors envision 
^1 sr sky coverage at ^mCrab sensitivity. For specificity, 
we adopt a FoV of 0.6 sr, a sensitivity of 1.5 mCrab in 100 s, 
and a bandpass 0.3 — 6 keV. Using the peak flux ^10~^ erg 
cm^^ s^^, the event could be detected out to z ~ 4. Based 
on a BAT detection rate of 1/7 yr^^, the WFXM rate would 
be ~1 — 2 yr~^. Currently operating X-ray /ASM instruments 
such as the RXTE/ ASM which has ^50 mCrab sensitivity in 
^6 X 10'' s (for blind searches) in the bandpass 2—12 keV 
would yield much lower rates, therefore their nondetections 
of these type of events are not unexpected. 

2.3. Falsifiability of the TDE interpretation 

In the picture of accretion by shredded stellar remnants, 
L cx t^^/'^ (Rees 1988), whereas in the accretion disk de- 



scription L cx t^^/^ for a super-Eddington disk with h/r ~ 1 
(CG09). The X-ray flux to date has shown a definite decrease, 
accompanied by large fluctuations. There are indications from 
comparing the current decay to the early light curve that a 
transition may have occurred around i ~ 6 d, which we ten- 
tatively interpret as the transition to disk decay. Therefore 
we should now (as of this writing, at i ~ TO + 100 d) be in 
the regime for which d\ogF/d\ogt ~ —4/3. The data to 
date do not definitely substantiate this prediction, but are gen- 
erally consistent. There are also large variations in the flux 
which may reflect jet instabilities (McKinney & Blandford 
2009) rather than variations in accretion rate onto the cen- 
tral engine. The present level of X-ray flux is about a factor 
10^ above the SwiftlXKY detection limit, therefore one could 
envision a more definitive decay law becoming manifest even- 
tually. If no such law emerges or if a different law emerges, 
that would argue against the TDE interpretation. If on the 
other hand the flux drops suddenly to an unobservable level, 
that may be consistent with a sudden change in the properties 
of the jet as the accretion rate drops below Eddington. Such a 
change is seen in the radio properties of X-ray binaries as they 
undergo state changes, e.g., from the high/thermal state to the 
low/hard state. For completeness, we note that a third option 
could be that nonlinear outcome of the TDE is completely dif- 
ferent from what we describe, in which case the nondetection 
of the t^^/'^ law would not be an argument against the TDE. 
In addition, even though the X-ray flux is currently dominated 
by the jet, it may be feasible to observe a thermal soft X- 
ray component from the disk. Using Equation (19), our disk 
model yields at i = 6 d a temperature at the ISCO (innermost 

stable circular orbit) GGMmbh/c^ of - 0.09 keV m^//^ 
which would produce a thermal component peaking at ^0.2 
keV and contributing ^10^^ in flux to the overall spectrum. 
According to our estimates a ^0.1 keV thermal component 
contributing only ^1% to the total observed emission could 
not be detected in the XRT spectra, which are dominated by 
the jet emission (Burrows et al. 2011). 

3. DISCUSSION AND CONCLUSION 

We have presented simple scalings laws for the post-TDE 
evolution. Three independent lines of reasoning point to a 
deeply plunging orbit i?p <C Rt if the TDE model is to apply 
to Sw 1644+57. 

(1) An application of the star fallback rate t^^^'^ to the 
early Sw 1644+57 light curve (as in Levan et al. 2011b) 
gives tfaiiback ^ 1 d. Thls implies that the orbit of the dis- 
rupted star must have been deeply plunging, ^ ~ 0.1. If 
Rp ss GA/mbh/c^ the commonly used Newtonian estimate 
for tfaiiback would need correcting (Ayal et al. 2000) because a 
standard Schwarzschild BH could not accommodate ^ < 0.2, 
for TOb7 = 1. This implies a maximal Kerr BH. Two first order 
effects would be Lense-Thirring precession and the general 
relativistic advance of the lines of apsides of ejecta trajecto- 
ries (Bardeen, Press, & Teukolsky 1972, Bardeen & Petterson 
1975). For 771^7 = 0.1 the minimum allowable ^ value (from 
Equation [3]) would decrease to ^0.04 for a Schwarzschild 
BH and ^0.02 for a Kerr BH, lessening the stringency im- 
posed by a ^ ~ 0.1 constraint. 

(2) The long term X-ray luminosity is ~'10^ — 10'^ times 
greater than expected for Eddington accretion. This points 
strongly to beaming, as has been inferred for GRBs. In ad- 
dition, in the long term fight curve there appears to be a 
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change in the decay properties at f ~ 6 d which we inter- 
pret as the onset of disk-dominated accretion. By comparing 
our theoretical decay rates for stellar debris fallback and ac- 
cretion in a freely expanding super-Eddington disk, we find 
that a small value ^ = Rp/Rx is required, 0.1 — 0.2, for 
-^^MBH — lO*" — 1O^M0. The observed X-ray luminosity at 
~2.5d Lx/^Eddington ~102 (lO^) for mm = 1 (0.1) would 
imply a beaming factor e;, ~ 10~^ (for f = ej = 0.1), similar 
to blazars (Sikora et al. 2005, Bottcher et al. 2007). The sud- 
den introduction of ^O.SMq of gas inside the ergosphere of a 
Kerr BH would launch powerful jets within a few local orbital 
time scales (McKinney 2005). Figure 1 of McKinney shows 
the steep dependence of jet efficiency on BH spin to mass ra- 
tio J/M, and in particular a steep upturn close to J/M ~ 1. 
In fact it may be that our putative constraint on Rp, namely 
Rp « GMmbh/c^ (which requires J/M ~ 1), is a necessary 
condition for this rare event. 

(3) The standard TDE rate 3 x 10~^ yr~^ galaxy"^ implies 
a Swift/BKT detection rate ~30 yr~^, for a beaming factor 
^10~^, which can be strongly excluded. A deeply plunging 
orbit, ^ ~ 0.1, lowers the TDE rate by a factor ^10^, giv- 
ing a rate in line with the BAT detection rate of ^1/7 yr^^. 
Our revised rate ^3 x 10^^ yr~^ galaxy^^ may need ad- 
justing to correct for relativistic effects. Ayal et al. (2000) 
compare Newtonian and relativistic TDE methods and find 
slight differences in i?p between the two sets of calculations, 
for the same initial trajectory. They also find that ^75% of 
the star becomes unbound, rather than ~ 50%. However, 
this depends strongly on the treatment of cooling, which was 
not modeled in their SPH (smoothed particle hydrodynam- 
ics) simulations. Another potential source of concern in our 
estimates is our use of the SS73-like algebraic scalings for 
the accretion disk. The SS73 methodology was expanded 
by Novikov & Thorne (1973) to include general-relativistic 
corrections including the full Kerr metric. Their multiplica- 
tive correction factors are based on Taylor expansions of the 
form 1 + ax + hx^ + cx^ + ... + a y + b + c where 
X = GMMBn/iRc^), y — J/M, and the factors a, b, c, 
a , b , c , ... are of order unity. These extra multiplicative 
terms would therefore modify our simple estimates by factors 
of order unity, even at small radii. 

It is worth noting that nearly all previous theoretical stud- 
ies of tidal disruption have, for simplicity, considered orbits 
for which the periastron radius Rp equals or is slightly inte- 
rior to Rt (Rees 1988, Evans & Kochanek 1989, CLG90, Ul- 
mer 1999, SQ09). Few works have considered deeply plung- 
ing disruptions (e.g., Guillochon et al. 2009). It appears Sw 
1644H-57 may require Rp <C Rt, and would be different in 
at least two important ways from the standard TDE. (1) First, 
the rates would be lower. The rate for deeply plunging or- 
bits would be less than that for traditional TDEs in which 
Rp ~ i?T by - {Rp/RtY, thus giving ~3 x lO"'^ yx'^ 
galaxy^ ^, if the more common normal, i.e., "total" TDE rate 
integrated over all Rp values is ^3 x 10^'"' yr^^ galaxy^^. 
(2) Second, to observe the event we need to be situated within 
the beam, which would only cover ^10^^ of the sky. Our 
estimate of the beam size is based on the ratio of observed X- 
ray luminosity to expected mass-energy accretion luminosity. 
There have been no detailed theoretical studies investigating 
the scenario we outline in this work. 

For an orbit in which Rp/Rp = 0.15, for example, 
the tidal force acting on the star at periastron would be 
^(0.15)"^ ~ 300 times that minimally needed to imbind 



the star The subsequent evolution might be more properly 
called "obliteration" rather than "disruption", hence we pro- 
pose the acronym TOE (tidal obliteration event). From SQ09 
and earlier studies, in the moving frame of the center-of-mass 
of the star, the acceleration of the outer layers due to the MBH 
tidal force ~(G'Mmbh/^p)(-R*/-Rp) acting over a dynam- 
ical time ~ {GMmbh/Rp)~^^^ results in velocity per- 
turbations to the stellar envelope Avp ~ v^{Rt/Rp)^^^, 
where is the escape speed from the stellar surface. For 
(Rp/Rp) — 0.15, the velocity perturbation exceeds the stel- 
lar escape speed by a factor ^20. The star would likely be 
completely shredded in ~10^ s. The half of the star rem- 
nants with orbits on inward trajectories would immediately 
establish a super-Eddington disk within the ergosphere of the 
MBH. 

The discovery of a new class of objects raises several in- 
teresting questions. For instance, why should a sudden accre- 
tion event onto a ~10^ — 1O^M0 BH resemble in any way 
a GRB? If there were an accretion event in which ^0.5Mq 
were suddenly introduced near the event horizon, the time 
scales would be consistent with what is observed. Within 
the first five days the light curve of Sw 1644H-57 showed 
flares with rise and decay times of '^lOO — 200 s. It is un- 
clear whether these variations stem from jet instabilities or 
are associated with disk time scales. An orbital time scale at 
tile tidal radius 2t: /^{Rt) = 1.5 hr r%^m~^'^ (0.5 hr for 
mt7 ~ 0.1). An after-the-fact examination of prior BAT data 
revealed a slow brightening over several days, significant at 
the ^2 — 2.5(7 level. The source was already bright in X-rays 
by the time XRT began observations. It appears we require a 
deeply plunging orbit, ^ ~ 0.1, and the orbital time 2'Kr/Vff, 
at the corresponding Rp is ^300 s for to;,7 = 1 (^100 s for 
= 0.1). The small periastron radius and general relativis- 
tic corrections such as the advance of the line of apsides would 
enormously speed up the TDE evolution compared to what is 
shown for instance in Figure 3 of SQ09. This is evidenced 
in the apparent -^1 — 10 hr timescale for the faUback time. 
The process responsible for utiUzing the accretion power in 
the inner disk to tap the BH spin, presumably the magnetoro- 
tational instability (Balbus & Hawley 1998) operating inside 
the ergosphere (McKinney & Narayan 2007ab), would pro- 
ceed with an initial growth rate given roughly by the local 
orbital frequency. Given a small seed field within the gas, 
it is likely that a large number >30 of e— folding times are 
required before the shear amplification of the local magnetic 
field leads to the nonlinear development of a strong jet. Weak 
precursor events may have contributed to low-level BAT ac- 
tivity preceding the main Sw 1644-h57 trigger. Considering 
the commonalities among jets observed to date in such a wide 
variety of objects. Nature appears to be telling us that the jet 
physics is largely independent of the mass of the central black 
hole. 
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